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ABSTRACT 

Context. The rapid neutron-capture process, or r-process, is known to be fundamental for explaining the origin of 
approximately half of the ^ > 60 stable nuclei observed in nature. In recent years nuclear astrophysicists have developed 
more and more sophisticated r-process models, by adding new astrophysical or nuclear physics ingredients to explain 
the solar system composition in a satisfactory way. Despite these efforts, the astrophysical site of the r-process remains 
unidentified. 

Aims. The composition of the neutron star outer crust material is investigated after the decompression that follows its 
possible ejection. 

Methods. The composition of the outer crust of a neutron star is estimated before and after decompression. Two 
different possible initial conditions are considered, namely am idealized crust composed of cold catalyzed matter and a 
crust initially in nuclear statistical equilibrium at temperatures around lO'^" K. 

Results. We show that in this second case before decompression and at temperatures typically corresponding to 8 10^ K, 
the Coulomb effect owing to the high densities in the crust leads to an overall composition of the outer crust in neutron- 
rich nuclei with a mass distribution close to the solar system r-abundance distribution. Such distributions differ, however, 
from the solar one due to a systematic shift in the second peak to lower values. After decompression, the capture of the 
few neutrons per seed nucleus available in the hot outer crust leads to a final distribution of stable neutron-rich nuclei 
with a mass distribution of 80 < ^ < 140 nuclei in excellent agreement with the solar distribution, provided the outer 
crust is initially at temperatures around 8 10^ K and all layers of the outer crust are ejected. 

Conclusions. The decompression of the neutron star matter from the outer crust provides suitable conditions for a 
robust r-processing of the light species, i.e., r- nuclei with A < 140. The final composition should carry the imprint of 
the temperature at which the nuclear statistical equilibrium is frozen prior to the ejection. 
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1. Introduction 

The r-process, or the rapid neutron-capture process, of stel- 
lar nucleosynthesis is invoked to explain the production of 
the stable (and some long-lived radioactive) neutron-rich 
nuclides that are heavier than iron and observed in stars of 
various metallicities, as well as in the solar system (for a 
review, sec Arnould et al. 2007). In recent years nuclear as- 
trophysicists have developed more and more sophisticated 
r-proccss models, trying to explain the solar system com- 
position in a satisfactory way by adding new astrophysical 
or nuclear physics ingredients. The r-process remains the 
most complex nucleosynthetic process for modelling from 
the astrophysics as well as nuclear-physics points of view. 
The site(s) of the r-process has (have) not been identified 
yet, since all the proposed scenarios face serious problems. 
Complex — and often exotic — sites have been considered in 
the hope of identifying astrophysical conditions in which 
the production of neutrons is high enough to give rise to a 
successful r-process. 

Progress in the modelling of type-II supernovae and 7- 
ray bursts has raised a lot of excitement about the so-called 
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neutrino-driven wind environment. However, until now a 
successful r-process cannot be obtained ab initio without 
tuning the relevant parameters (neutron excess, entropy, 
expansion timescale) in a way that is not supported by 
the most sophisticated existing models. Although these sce- 
narios remain promising, especially in view of their poten- 
tial to significantly contribute to the galactic enrichment 
( Argast et al. 2004) , they remain handicapped by large un- 
certainties associated mainly with the still incompletely un- 
derstood mechanism that is responsible for the supernova 
explosion and the persistent difficulties obtaining suitable 
r-process conditions in self-consistent dynamical explosion 
and neutron-star cooling models (Hiidepohl et al. 2010; 
Fischer et al. 2010). In addition, predictions of the detailed 
composition of the ejected matter remain difficult owing to 
the remarkable sensitivity of r-process nucleosynthesis to 
uncertainties of the ejecta properties. 

Early in the development of the theory of nucleosynthe- 
sis, an alternative to the r-process in high-temperature su- 
pernova environments was proposed (Tsuruta et al. 1965). 
It relies on the tendency of matter at high densities (typ- 
ically p > lO^'^ gcm^'^) to be composed of nuclei lying on 
the neutron-rich side of the valley of nuclear stability as a 
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result of endothermic free-electron captures. This so-called 
'neutronization' of matter is possible even at zero temper- 
ature. The astrophysical plausibility of this scenario in ac- 
counting for the production of the r-nuclides has long been 
questioned. It remained largely unexplored until the study 
of the decompression of cold neutronized matter resulting 
from tidal effects of a black hole on a neutron star (NS) 
companion (Lattimer et al. 1977; Meyer 1989). Recently, 
special attention has been paid to NS mergers following 
the confirmation by hydrodynamic simulations that a non- 
negligible amount of matter can be ejected (Janka et al. 
1999; Rosswog et al. 2004; Oechshn et al. 2007). The ejec- 
tion of initially cold, decompressed NS matter might also 
happen in other astrophysical scenarios like giant flares in 
soft-gamma repeaters, the explosion of an NS eroded be- 
low its minimum mass (e.g. Sumiyoshi et al. 1998), or the 
equatorial shedding of material from very rapidly rotating 
supramassive or ultramassive NSs (see Arnould et al. 2007 
for more details). 

The composition of the NS matter of the inner crust, be- 
fore and after decompression, has already been studied in 
detail in Goriely et al. (2005) and Arnould et al. (2007). It 
was found that the final composition of the material ejected 
from the inner crust is expected to depend on the initial 
density, at least for the upper part of the inner crust at 
Pdrip < /o[g/cm3] < 10^2 (where pdrip - 4.2 lO^g/cm^ is 
the neutron drip density) . For the deeper inner-crust layers 
(p > 10^^ g/cm^), large neutron-to-seed ratios drive the nu- 
clear flow into the very heavy mass region, leading to flssion 
recycling. As a consequence, the abundance distribution is 
now independent of the initial conditions, especially the ini- 
tial density. In both cases, the abundance distribution was 
found to be in close agreement with the solar distribution 
for A > 130 nuclei (Goriely et al. 2005; Arnould et al. 2007). 

Although the outer crust is far less massive than the 
inner crust, the ejection of the inner crust cannot take 
place without at the same time leading to the ejection of 
at least some outer crust material. The outer crust typi- 
cally amounts to 10~^ to 10"'' Mq, depending on the NS 
mass and radius (Pearson et al. 2011) so that depending on 
the frequency and fraction of its possible ejection, it may or 
may not contribute significantly to the galactic enrichment. 
We will not get into these considerations in the present pa- 
per, but restrict ourselves to estimating the composition of 
the outer crust material before (Sect. [2]) as well as after 
ejection (Sect.[3|, assuming that nuclear statistical equilib- 
rium (NSE) could be established before ejection occurred. 
We shall consider two cases: the first one corresponds to an 
initially cold NS crust and the second one to an NS crust 
initially in NSE at a temperature of the order of 10^° K at 
which nuclear statistical equilibrium (NSE) is responsible 
for setting the composition prior to the ejection. 



2. Composition of the NS outer crust in equilibrium 

The composition of the outer crust needs to be estimated 
in each layer, each one being characterized by a given den- 
sity and pressure. The overall abundance distribution for 
the outer crust (0 < p{r) < Pdnp) can then be integrated 
over a pressure column P{r) (where r is the radial coor- 
dinate) given by the Tolman-Oppenheimer-Volkoff (TOV) 



equations (considering non-rotating NS) 

dP _ {M{r) + 4Trr^P{r)/c^}{£{r) + P(r)} 

dr ^ r{rc^ - 2GMir)} ^ ' 

with 

^ = 47rf (r)rVc2 , (2) 
dr 

in which £ denotes the total energy density, including rest 
mass (in the outer crust, £ ~ pc^) and Ai is the gravita- 
tional mass. We choose here a standard NS with a mass of 
1.5 Moand 13 km radius. Since the equation of state (EoS) 
and the composition remain sensitive to the nuclear physics 
ingredients adopted (Pearson et al. 2011), we will consider 
here three different Hartree-Fock-Bogolyubov (HFB) mass 
models, namely the Skyrme-HFB mass models HFB-19 and 
HFB-21 (Goriely et al. 2010) and the Gogny-HFB mass 
model DIM (Goriely et al. 2009). For each mass model the 
EoS has been consistently taken into account in solving the 
TOV equation. In the three cases, a baryonic mass of about 
5 10~^ M0 (Pearson et al. 2011) is found in the outer crust 
of such an NS and about 90% of the outer-crust mass is lo- 
cated in the density range of 6 10^" g/cm^ ^ P l£ Pdrip- For 
densities above the drip density, free neutrons are present 
and their contributions to the EoS should be included. A 
transitional regime into the inner crust can be included 
in the calculation by estimating the pressure and free en- 
ergy associated with the corresponding gas of free neu- 
trons obtained self-consistently with the effective interac- 
tions BSkl9, BSk21 or DIM corresponding to the above- 
mentioned mass tables. 

2.1. The cold NS crust 

A fairly common picture of an NS (Baym et al. 1971; 
Pethick & Ravenhall 1995) is that they consist of "cold cat- 
alyzed matter", i.e., electrically neutral matter in its abso- 
lute ground state in the sense of complete nuclear and beta 
equilibrium at temperature T — 0. Detailed calculations 
of the composition of the NS material initially at a den- 
sity below the drip density {pdnp — 4.2 x lO"'^ gcm~^) have 
been performed over the years (I3aym et al. 1971; Pethick & 
RavenhaU 1995; Haensel & Pichon 1994; Riister et al. 2006; 
Pearson et al. 2011). We have repeated the calculation of 
Baym et al. (1971) with the updated nuclear physics inputs 
to determine the outer-crust composition. This calculation 
minimizes the free Gibbs energy per nucleon at T = 0. All 
details can be found in Pearson et al. (2011) and will not 
be repeated here. 

For densities above 3 lO^g/cm^, the outer crust is es- 
sentially made of iV = 50 and iV = 82 neutron-rich nu- 
clei. More precisely, for 10^ < p[g/cm^] <, 5 10^*^, we find 
TV = 50 nuclei with 80 < A < 86. At these densities, only 
nuclei with experimentally known masses are involved. This 
is not the case at higher densities, i.e., in the most mas- 
sive part of the outer crust, where we use the HFB-19, 
HFB-21 or DIM mass tables to complement experimental 
masses. For 10^^ < p[g/cm'^] < pdrip, iV = 82 nuclei with 
120 < A < 126 populate the outer crust. 

The variation of the T = composition of the outer 
crust as a function of density is illustrated in Fig. [T] In 
this density range, the electron fraction Ye varies from 0.43 
down to 0.31 at the drip density. Once in possession of the 
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Fig. 1. Ground-state composition (charge, neutron and 
mass numbers as well as free neutrons) of the outer crust 
and of the shallow layers of the inner crust as a func- 
tion of the density. Predictions with HFB-19 masses (solid 
lines) (Goriely et al. 2010) are compared with those ob- 
tained with the DIM masses (dotted lines) (Goriely et al. 
2009). Experimental masses (Audi et al. 2003, 2010) are 
used whenever available. 



solution p{r) to the TOV equations, it is easy to calcu- 
late the baryonic mass of any shell whose inner and outer 
radii ri and r2, respectively, are chosen to correspond to 
given densities p{ri) and p{r2) (Pearson et al. 2011). The 
abundances of the different nuclides in the outer crust can 
then be read off from Fig. [l] the results being shown in 
Fig. [2] Slightly different distribut ions are obtained with 
different mass models, but generally, because of the high 
mass included within the high-density region (close to the 
drip point), nuclides with A ~ 80 or 118 < A < 126 will 
be dominant, as can be seen in Fig. [2j 

If a decompression of the T — outer crust mate- 
rial takes place, only /3-decays towards the valley of sta- 
bility are expected, no free neutrons being present (ex- 
cept through the possible /3-delayed neutron emission). 
Consequently, after decompression only a restricted distri- 
bution of 78 ^ A ^ 86 and 120 ^ A ^ 126 nuclei can be 
expected, as discussed in Sect. [3T| 

2.2. The hot NS crust 

The T — picture of the NS crust is clearly an idealization, 
as it implies that the star has had an infinite time to cool 
down and maintain (or restore) thermodynamic equilib- 
rium since its creation in the aftermath of a gravitational- 
collapse supernova explosion. In reality, not only will the 
actual temperature of any particular NS be non-zero, but 
the equilibrium configuration corresponding to a still higher 
temperature might have become "frozen in" . The outer NS 
crust may in fact have a very different composition at non- 
zero temperature due to the specific softness of the distri- 
bution of the Gibbs free energy per nucleon. 

As shown in Fig. [3] many nuclei yield a free energy per 
nucleon that differs from the equilibrium one by no more 
than a few tens of keV. At finite temperatures, thermal 
fluctuations could therefore significantly broaden the dis- 
tribution of equilibrium nuclides. At Tg <^ 4 — 5 (where Tg 



Fig. 2. Abundance distribution integrated over the whole 
outer crust for an initially cold NS of mass 1.5 M0and 
13 km radius. The distributions are arbitrarily normalized 
to the solar system r-abundance distribution, shown for 
comparison. The distributions are obtained with the HFB- 
19 (squares), IIFB-21 (circles) and DIM (diamonds) mass 
models when experimental masses (Audi et al. 2003, 2010) 
are not available. 
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Fig. 3. Distribution of the Gibbs free energy per nucleon 
at T = and p = 3.7 10^^ g/cm^. AG corresponds to the 
difference of the free energy per nucleon with respect to the 
minimum value obtained at this density for Z = 38, N = 88 
and A = 124. The distribution is obtained with the IIFB-19 
mass model (Goriely et al. 2010) when experimental masses 
are not available. The open squares represent stable nuclei. 



is the temperature expressed in lO^K), an (n,7)-(7,n) equi- 
librium may be established within the isotopic chain corre- 
sponding to the most probable Z value at the given density. 
For higher temperatures, an NSE could be reached and the 
abundance of a given nucleus i oi Z protons and N neu- 
trons {A = Z + N) is now given by (Bravo & Garcia-Senz 
1999; Nadyoshin & Yudin 2005; Arcones et al. 2010) 



Y,, = 



uji ( kTArriu 



3/2 



(3) 
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where Ui is the T-dependent partition function, and 

Qi = [Zmp + Nm,i - mi] <? + Z[ic,p - Mc(2', A) , (4) 

includes the binding energy, expressed in term of the nu- 
cleon and nuclear masses Mi , as well as the Coulomb correc- 
tions to the chemical potential [ic arising from the Coulomb 
contribution to the free energy, which becomes significant 
for heavier nuclei. Eq.|3]assumes that nuclei follow Maxwell- 
Boltzmann statistics, but the presence of the degeneracy 
parameters 7y„ and ryp takes into account the possibility 
that nucleons follow the more general Fermi-Dirac statis- 
tics, which will certainly be the case for free neutrons close 
to the drip density. These parameters are related to the 
number density Uq through the usual expression 



87rV2 



(5) 



where /3 = kT/mc^ is the relativistic parameter and are 
the Fermi functions 



dx . 



(6) 



More details on the high-density NSE equations can be 
found in Bravo & Garcia-Senz (1999) and Arcones et al. 
(2010) 

To estimate the composition of the hot crust, the same 
EoS is used as in Pearson et al. (2011), though the T- 
dependent non-ideal Coulomb interaction between ions as 
well as between ions and electrons and between electrons in 
the approximation of a rigid electron background is taken 
from Haensel et al. (2007) and DeWitt et al. (1996). In this 
approximation, the Coulomb chemical potential of species 
i is given by 



(7) 



where fc is the Coulomb free energy per ion in units of 
kT. We would like to point out that we did not include 
here the contribution of fcT/3 x dfc/d\\iTi as proposed by 
Glazyrin & Blinnikov (2010), since this contribution implies 
a variation of the electron density and should therefore not 
be included in the NSE equations which are calculated at 
a fixed value of the electron density (see the discussion p. 4 
of Glazyrin & Blinnikov 2010). 

For a Coulomb liquid, fc can be expressed as (Haensel 
et al. 2007) 



c 



riikT 



= Aiv/r.(A2 + r,) 



(8) 



-Al X A2 In (yfyl^ + v/TTrv^) 



-2A. 



arctan 



F, - B2 In 1 



B. 



In 1 



Bi 



with Al = -0.9070, A2 = 0.62954, ylg 0.27710, Bi = 
0.00456, B2 = 211.6, S3 = -0.0001 and B4 = 0.00462. The 
corresponding contribution to the internal energy and the 
pressure can be found in Haensel et al. (2007). 

The Coulomb liquid approximation is adopted since at 
the temperatures considered here (typically Tg = 5 — 10), 
the Coulomb coupling parameter, 
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Fig. 4. Melting temperature T™ as a function of the density 
p, estimated from Eq. (ToJ) and the T — composition 
shown in Fig.[l]using HFB-19 masses when no experimental 
masses are available. 



where is the ion-sphere radius, remains smaller than the 
melting value F,„ = 175.0±0.4 (Pothekin & Chabrier 2000). 
Equivalently, the temperature T is higher than the melting 
temperature. 
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1/3 



K, 



(10) 



as shown in Fig. |4j where Tm has been evaluated on the 
basis of the T — composition given in Fig. [l] (based on 
HFB-19 masses). 

To estimate the NSE composition of a hot NS crust 
still requires the determination of the distribution of the 
electron fraction Ye, i.e., its variation with the density n. 
We assume here that the NS cooled sufficiently slowly for 
/3-equilibrium to be maintained down to a temperature of 
about T ~ 10^" K, at which point the Ye distribution be- 
comes frozen in as the NS continues to cool to lower tem- 
peratures. However, we assume that NSE is maintained 
down to t empe ratures of the order of T ~ 7—10 10^ K 
(see Sect. 3.2 for a discussion of the NSE timescales). 
The Ye distribution corresponding to /3-equilibrium is ob- 
tained by minimizing the free Gibbs energy per nucleon 
at r = 10^° K; for the physics used, see also Pearson et 
al. (2011). We use the traditional linear mixing rule (also 
known as the additive approximation) to estimate the total 
free Gibbs energy of the multicomponent plasma by weight- 
ing the individual free energy of each species by its molar 
fraction in the mixture. With respect to the T = case, the 
Ye distribution is found to be slightly higher, and obviously 
much smoother, as shown in Fig. [5j The implied decrease 
in the neutron fraction occurs despite the liberation of free 
neutrons at non-zero temperatures. 

The NSE abundance distributions for the whole outer 
crust at four different temperatures (between 7 and 10 GK) 
and integrated over a pressure column given by the TOV 
equations are shown in Fig. [6j The outer NS crust in NSE 
at a temperature of Tg ~ 7—10 is seen to be made of 
r-nuclei with a distribution similar to that of the solar sys- 
tem, although the second peak is slightly shifted to lower 
masses around A ~ 126, as is also found for the T = 
composition (Fig. [2| . The second peak is also found to be 



S. Goriely et al.: The r-process nucleosynthesis in neutron stars 



5 



0.45 > 



0.40 



0.35 



0.30 



0.25 




10" 



10' 



U 10-2 

-a 

3 

< 

> 

Pi 



10 



10- 



10' 



= 1 1 1 1 1 1 

; Solar 


1 1 1 1 1 1 1 1 1 


1 ' ' ' 1 ; 

-■-HFB-19 ; 






-^HFB-21 1 






— ^DIM 














1 


1 

1 . . . 1 . . . 





80 



P [MeV/fm ] 



Fig. 5. Ye distribution in the outer crust, assuming 
(neutrino-less) /3-equilibrium, as a function of the pressure. 
The dashed line corresponds to T = and the solid line to 
T = 101° K. 



Fig. 7. Same as Fig. [b] for Tg = 8 and different mass mod- 
els, namely HFB-19, HFB-21 (Goriely et al. 2010) or DIM 
(Goriely et al. 2009). 
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Fig. 6. NSE abundance distribution integrated over the 
whole outer crust for 4 different temperatures assuming 
the initial Yg distribution is determined by /3-equilibrium at 
T = 10^° K. The distribution is obtained with the HFB-19 
mass model (Goriely et al. 2010) when experimental masses 
are not available. The solar system r-abundance distribu- 
tion (dotted circles) is shown for comparison. 



rather pronounced with respect to the first peak for low 
temperatures (Tg = 7), but to be smoothed away, on the 
other hand, for high temperatures (Tg = 10). The abun- 
dance distribution is sensitive to the adopted mass models, 
as illustrated in Fig. [7] but is otherwise extremely robust, 
since as far as the nuclear physics is concerned it depends 
only on binding energies. The Coulomb correction plays a 
vital role by shifting significantly the abundance distribu- 
tion towards the high-mass region. In particular, no ~ 82 
would be found in the crust if the Coulomb interaction was 
omitted in the NSE equations, as shown in Fig. |8] 

3. Decompression of the NS outer crust 

Without focussing clearly on a particular ejection scenario 
of outer-crust material of an NS, we decribe the density evo- 
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Fig. 8. Same as Fig.[7]including and excluding the Coulomb 
interaction (the HFB-19 mass model is used). 



lution of an ejecta clump in simple, but as general terms as 
possible. For setting the values of the free parameters we 
will refer to existing results from simulations of the dynam- 
ical merging of binary NSs or an NS with a black hole. In 
these events near-surface matter can be tidally stripped and 
ejected, and we take the ejection dynamics of such events 
as template for our expansion model. 

The evolution of the matter density is approximated 
by a one-zone model of the pressure-driven expansion of 
a self-gravitating clump under the influence of tidal forces 
(exerted by the compact source object of the ejecta) along 
the escape trajectory of the clump. In order to determine 
the time evolution of the clump radius R{t), we integrate 
numerically the Newtonian equation of motion 



d^i? 
dt^ 



R 



(11) 



Here the first term corresponds to the confining acceler- 
ation by the self-gravity of the clump, approximated by 
-GMJR^ w -^GpR when w ^pR^ is the clump 
mass with average density p. The inflating effect of the gas 
pressure P of the clump is repesented by the second term 
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in Eq. (11), where we use a one-zone representation of the 
pressure gradient according to /3~^dP/dr « P/{pR). The 
third term describes the acceleration corresponding to the 
tidal stretching by the external gravitational force of the 
NS, 



GAL 



tidal 



{rc + Rf 



2GMnsR 



(12) 



where M^s is the NS mass and rc{t) the time-dependent 
radial position of the clump center. In the transforma- 
tion leading to the last term in Eq. (12 1 we made use of 
the relation R <^ rc- Assuming that the ejection velocity 
Vc = drc/dt of the clump is given by the escape velocity on a 
parabolic orbit, i.e., = GMns/rc, we can solve the equa- 
tion of motion of the clump center, dvc/dt — y^2GAIns/rc, 
to obtain 



3 t 

2 Tnoi 



2/3 



(13) 



with To being the radial position of the clump center at 
the beginning of the ejection, t^sc = ro/vcsc the escape 
timescale, and Wcsc = \/2GMns/ro the escape velocity. 
Using rc(i) of Eq. (Il3| in Eq. (pi) we get 



d^R 
~diP~ 



R 



tidal 



(14) 



which appears as the last term on the rhs of Eq. (111. 
Analytic estimates for the conditions at the surface of an 
NS in agreement with numerical results of NS merger sim- 
ulations (Ruffert & Janka 2001) yield values of the escape 
timescale Tgsc in the range of 10~^ to 3 10^* s. The value 
of Tesc — 3 10~^ s is adopted in this work. 



In solving the simple model equation, Eq. (11), the ini- 
tial clump radius Rq plays an important role, because it de- 
termines the initial acceleration of the clump expansion and 
therefore its expansion timescale. It can be used as a free 
parameter, leading to different density decline rates during 
decompression. Here we choose a value of i?o ~ 2 km, which 
is found to reproduce the density evolution of the major- 
ity of NS merger ejecta (Ruffert & Janka 2001) fairly well. 
The pressure P of the clump medium is determined from 
the EoS of Timmes & Arnett (1999) and is consistently 
evolved for the changing conditions of density and temper- 
ature (including possible /3-decay heating) during the de- 



compression history. The integration of Eq. (11) yields the 
time-dependent clump radius R{t), with which (for con- 
stant clump mass) the clump density can be computed as 
p{t) = po{Ro/R{t)f. 

The composition change during the decompression is 
followed with a full network, where in particular the (3- 
decay processes have been taken from the updated version 
of the Gross Theory (Tachibana et al. 1990) and neutron 
capture rates are consistently estimated with the TALYS 
code (Goriely et al. 2008) on the basis of the nuclear mass 
model used for the initial conditions. The temperature evo- 
lution is followed as described in Meyer (1989) and Goriely 
et al. (2005) on the basis of the laws of thermodynamics and 
possible nuclear heating through /3-decay processes. The de- 
compression along the above-defined trajectory of both the 
cold NS outer crust and the outer crust initially in NSE are 
studied below. 
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Fig. 9. Same as Fig. [i] after decompression. 



3.1. The cold NS crust 

The final isobaric abundance distribution of the matter 
ejected from the cold NS crust initially assumed to be at 
/3-equilibrium can be expected to remain relatively similar 
to the one prior to the ejection, since only /3-decay and 
/3-delayed neutron emission (with the possibility of recap- 
turing the emitted neutrons) may change the initial compo- 
sition. The final isobaric abundance distribution is shown 
in Fig.[9j the initial composition being given in Fig. [2]prior 
to the decompression. Globally, the final distribution is ob- 
viously far from matching the overall solar system pattern. 
Some differences are seen depending on the nuclear physics 
ingredients adopted. Interestingly, the decompression of the 
cold NS outer crust appears as a potential site for produc- 
ing almost exclusively the 115 <, A <, 124 r-nuclei. In some 
site-independent parametrized r-process models, these nu- 
clei have been found to be underproduced, though this con- 
clusion is obviously subject to large nuclear physics and 
astrophysics uncertainties (for a review, see Arnould et al. 
2007). 

3.2. The hot NS crust 

As discussed in Sect. |2.2[ the outer crust initially in NSE 
at temperatures around Tg = 7—10 differs considerably 
from the cold T = case. In addition, for layers close to 
the drip density, free neutrons are present in significant 
numbers due to the high initial temperatures. In this case, 
a few neutrons per seed nucleus are available and may be 
captured during the expansion phase so that the initial NSE 
mass distribution is modified. The abundance distributions 
in a layer with initial pressure Pq = 3 10^** MeV/fm^ and 
density po = 2.7 lO^g/cm^, initially in NSE at Tg = 8 is 
shown in Fig. [TOj The evolution of the temperature, density 
as well as the corresponding radioactive power due to /3- 
decays is shown in Fig. [TT] fo r the same ejected layer as 
the one considered in Fig. |10[ In particular, it can be seen 
that the /3-decay heating slows down the temperature drop 
already a few ms after ejection. 

After decompression, due to the capture of about 5 free 
neutrons per seed nucleus, the distribution is shifted to- 
wards the A ~ 130 peak and shaped by the lower tem- 
peratures found at the time of the neutron captures. The 
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A 

Fig. 10. Abundance distribution before and after de- 
compression of a layer with initial pressure Pq = 
4 10"^ MeV/fm^ and density po = 3.4 lO^g/cm^ and ini- 
tially in NSE at Tg = 8. The calculation was performed 
with the HFB-19 masses and corresponding reaction rates. 
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Fig. 14. Time needed to reach NSE as a function of the 
temperature for different values of Ye at a constant density 
of p = 3 10^^ g/cm~'^. Symbols corresponding to the same 
value of Ye are connected by a solid line. Ye — 0.50 is shown 
by squares, 0.45 by circles, 0.40 by triangles and 0.33 by 
diamonds. 



10 




Fig. 11. Evolution of the temperature, density and ra- 
dioactive heating rate per unit mass {Q) resulting from 
the decompression of the layer with an initial pressure 
Pq = 4 10-'* MeV/fm^ and density po = 3.4 lO^g/cm^ 
and initially in NSE at Tq = 8 10^ K. The inftial radioac- 
tive decay heat Qq amounts to 6.6 10*® erg g"* s~^. 



final distribution is found to be independent of the details 
characterizing the expansion. In particular for a faster ex- 
pansion obtained with the value of Tosc = lO"'* s a virtually 
identical distribution is found. 

The abundance pattern for the outer crust integrated 
up to a maximum pressure P is given in Fig. [l2] for the 
initial temperature of Tg = 8. If the whole outer crust 
is considered {P < Pdrip = 5 10"^ MeV/fm'^), an over- 
all agreement with the solar distribution is obtained in the 
whole 80 <^ A ^ 150 mass region. An overproduction of the 
A > 140 r-nuclei (relative to the solar distribution) could 
be obtained if a sizable part of the inner crust was ejected 
at the same time (Arnould et al. 2007). When the mass 
ejection is restricted to parts of the outer crust (up to a 
pressure value P < Pdnp), different r-abundance distribu- 
tions are obtained, as shown in Fig. |12| The production of 



r-nuclei in such a scenario is clearly sensitive to the thick- 
ness of the outer crust that is ejected and therefore depends 
on the ejection mechanism that is invoked. 

The overall abundance distribution depends not only 
on the fraction of the ejected crust but also on the initial 
temperature at which the NSE has been frozen in. We show 
in Fig. [TSjthe abundance distributions obtained for different 
initial temperatures between 7 and 10 10^ K, assuming that 
the whole outer crust is ejected. 

It could be seen as purely fortuitous that temperatures 
around Tg ~ 8 give rise to r-abundance distributions in 
agreement with the solar distribution. However this is in 
agreement with the astrophysical scenario considered here. 
In particular, the temperatures considered here typically 
correspond to those at which NSE can be dynamically 
achieved in cooling events. We have calculated the typi- 
cal timescale tjsise needed to reach NSE for a density of 
p = 3 10** g/cm~^ (typical of the most contributing part 
of the outer crust) and for different values of Ye (as found in 
the outer crust). This timescale is obtained by estimating 
the time needed to reach a steady state at constant temper- 
ature and density considering different initial seed nuclei for 
the chosen value of Ye. A full network including all strong 
and electromagnetic interactions (but no weak interactions) 
is used for this purpose. As shown in Fig. [MJ for matter 
with initial values of Ye = 0.33 — 0.40, it takes about 1 to 
20 ms to reach an NSE at Tg = 8 and p = 3 10** g/cm'^, 
while at Tg ~ 9, around 0.2 ms are required for the most 
neutron-rich conditions (Ye = 0.33). Such ms timescales 
are characteristic of dynamical scenarios of interest here 
for the potential mass ejection (e.g. in the bursts of soft 
gamma-repeaters, during NS mergers, . . . ), so that prior to 
the ejection the NSE should be achieved at temperatures 
typically above Tg ~ 8 — 9 and can be expected to be frozen 
below during the ejection. The final production of r-nuclei is 
clearly sensitive to the temperatures at which NSE freezes 
in so that depending on the dyn amical timescales different 
distributions (as seen in Fig. 13 1 are found. 
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Fig. 12. Final abundance distributions for the whole outer crust (as integrated up to a pressure P given in the legend in 
MeV/fm'^) after decompression when initially in NSE at Tg = 8. The calculations are performed with the HFB-19 mass 
model and corresponding rates. The distributions are compared with the solar r-abundance distributions (dotted circles) . 
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Fig. 13. Final abundance distributions of the outer crust material after decompression, if initially in NSE at different 
temperatures ranging between Tg = 7 and Tg = 10. The calculations are performed with the HFB-19 mass model and 
corresponding rates. The distributions are compared with the solar r-abundance distributions (dotted circles). 
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Fig. 15. Same as Fig. 13 for an initial NSE temperature 
of Tg = 8 and five different mass models (and correspond- 
ing reaction rates), namely the HFB-19, HFB-21 and DIM 
mass models (upper panel) and the Duflo & Zuker (DZ; 
1995) and FRDM (MoUer et al. 1995) models (lower panel).' 



Uncertainties in the nuclear physics input can influ- 
ence the calculation of the final abundances. These con- 
cern essentially nuclear masses and the corresponding reac- 
tion rates. We show in Fig. [15] the abundance distributions 
obtained with NSE distributions and reaction rates deter- 
mined on the basis of different mass models. The initial NSE 
temperature is 8 10^ K and the whole outer crust mass 
is assumed to be ejected. As shown in Fig. [15] the exact 
strength, location, and width of the = 50 and = 82 r- 
process peaks are sensitive to the mass model adopted. The 
production of the 90 < yl < 120 r-nuclei can also be signifi- 
cantly modified according to the nuclear physics input. The 
importance of nuclear masses in this specific r-process nu- 
cleosynthesis scenario is essentially linked to the initial NSE 
conditions, which define not only the most abundant species 
initially present in the outer crust, but also the amount 
of free neutrons available for the neutron-capture process 
during the decompression. In contrast, the impact of un- 
certainties affecting the /3-decay rates is found to be rather 
limited (as illustrated in Fig. 16 1, though /3-delayed neutron 
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emission plays some role in smoothing the final abundance 
pattern. 



Fig. 16. Same as Fig. 13 for an initial NSE temperature 
of Tg = 8.5 and two different calculations of the /3-decay 
rates, namely the Gross Theory (version 2; Tachibana et 
al. 1990) and the Tamm-Dancoff approximation (TDA) of 
Klapdor et al. (1984). The calculations were performed with 
the HFB-21 mass model and corresponding rates. 



4. Conclusion 

Most of the problems faced in understanding the origin 
of r-process elements and observed r-abundances are re- 
lated to our ignorance of the astrophysical site that is ca- 
pable of providing the required high neutron flux. We have 
shown here that the decompression of the NS matter from 
the outer crust provides suitable conditions for a robust r- 
processing of the light species, i.e., r-nuclei with A < 140, 
provided the outer crust is initially in NSE at temperatures 
around Tg ~ 8 and the ejection mechanism allows for the 
full outer crust to be expelled to the interstellar medium. 
During the decompression, the free neutrons (initially lib- 
erated by the high temperatures) are re-captured leading 
to a final pattern similar to the solar system distribution. 

While the total mass included in the outer crust only 
amounts to some 10^^ — 10"'' Mq, its ejection leads to 
a galactic enrichment in stable neutron-rich nuclei with 
a final composition close to what is now observed in the 
solar system. The final composition should carry the im- 
print of the initial temperature at which the NSE is frozen 
in prior to the ejection, as well as the density region of 
the outer crust that gets ejected. The final abundances are 
also affected by nuclear uncertainties, most particularly nu- 
clear masses. The similarity between the predicted and solar 
abundance patterns as well as the robustness of the predic- 
tion against variations of input parameters (such as expan- 
sion timescales or initial Yg distribution) make this site one 
of the most promising that has been proposed, deserving 
further exploration with respect to various aspects such as 
nucleosynthesis, hydrodynamics, and galactic chemical evo- 
lution. 
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